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Abstract. We give a short survey on computational techniques which can 
be used to solve the representation conversion problem for polyhedra up to 
symmetries. We in particular discuss decomposition methods, which reduce 
the problem to a number of lower dimensional subproblems. These methods 
have been successfully used by different authors in special contexts. Moreover, 
we sketch an incremental method, which is a generalization of Fourier-Motzkin 
elimination, and we give some ideas how symmetry can be exploited using 
pivots. 



1. Introduction 

By the Farkas-Minkowski-Weyl Theorem a convex polyhedron in IR^ has two 
representations. It can either be described by a finite set of linear inequalities 
(facets) or by a finite set of generators (vertices and rays). Precise definitions are 
given in Section [2] 

One of the most fundamental problems in the theory of polyhedra and its ap- 
plications, such as Combinatorial Optimization or Computational Geometry, is the 
conversion between two different descriptions. Many algorithms for this representa- 
tion conversion have been proposed (see for example |MR80j . |ABS97] , [BFM98 , 
|JZ04j ). For certain classes of polyhedra efficient methods are known, but there is 
no approach known which efficiently solves the problem in general. Programs like 
cdd [ Fu95] . Irs |Av93j . pd [Ma97j . porta |ChLo97] (or polymake jGJOOj either 
relying on some of the others or using its own method closely related to cdd) allow 
conversion of the representation of a polyhedron. Since the programs are imple- 
mentations of quite different methods, their efficiency may vary tremendously on a 
given example. 

Many interesting polyhedra both pose difficulties for the standard representa- 
tion conversion approaches and have many symmetries that could potentially be 
exploited. In many applications it is sufficient (or at least necessary) to obtain a 
list of inequalities or generators up to symmetries. In the present paper we give a 
brief survey of approaches that can be used for the representation conversion prob- 
lem up to symmetries. We do not discuss their asymptotic complexity (which, in 
the worst case, is not encouraging), but rather refer to previous papers where the 
methods have proven themselves on difficult instances that could not otherwise be 
solved. For the new approach discussed in Section [7] we provide some experimental 
data ourselves. 

The paper is organized as follows. In Section [2] we give some basic notations 
and facts from the theory of convex polyhedra. In Section [3] we consider different 
notions of symmetries and describe how they can practically be obtained using a 
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graph automorphism computation. In Section [4] we describe the group theoretical 
notions used in the representation conversion methods discussed in the remainder 
of the paper. In Section [5] we consider decomposition methods which reduce the 
given problem to a number of smaller problems. These approaches have been used 
quite successfully by different authors. In Section [6] we describe the incremental 
Cascade algorithm and in Section [7] we show how symmetry can be exploited in a 
simplex pivot based algorithm. 

2. Convex polyhedra 

In this section, we give a brief introduction to some basic concepts and termi- 
nology of (convex) polyhedra. For more details on polyhedra, we refer to the books 
[Zi98] , |Gr03j . |Sch86j . 

Given the vector space M d , denote by (M. d )* its dual vector space, i.e. the vector 
space of linear functionals on M. d . A convex polyhedron V C M. d can be defined by a 
finite set of linear inequalities 

V = {x eR d : fi{x) > h,i = 1, . . . ,m} 

with fi G (M. d )* and bi G R for i — 1, . . . , m. If the number of inequalities m in the 
description is minimum, we speak of a non-redundant description. The dimension 
dim'P of V is the dimension of the smallest affine subspace containing it. Under 
the assumption that V is full-dimensional, i.e. dim? 3 = d, every inequality i of a 
non-redundant description defines a facet {x G V : fi(x) = bi} of V, which is a 
(d — l)-dimcnsional convex polyhedron contained in the boundary of V . 

By the Farkas-Minkowski-Weyl Theorem (see e.g. |Sch86] . Corollary 7.1a), V 
can also be described by a finite set of generators: 

V = conv{«i,...,i>fc}+cone{u fc+ i,...,w„} 

n k 

i=l i=l 

where Vi G K d for i — l,...,n. If the number of generators is minimum, the 
description is again called non-redundant. In the non-redundant case, the generators 
Vi, i = 1, . . . , k, are called vertices and M>oWi, i = k + 1, . . . , n, are the extreme rays 
of V . In case V is bounded we have n = k and we speak of a convex polytope. 

The representation conversion from a minimal set of generators into a minimal 
set of linear functionals (or vice versa) is called the dual description problem. By 
using homogeneous coordinates, the general inhomogeneous problem stated above 
can be reduced to the homogeneous one where bi — and k — 0. For example, we 
embed V £ M. d in the hyperplane Xd+i = 1 in M. d+1 and consider the closure of its 
conic hull. The so-obtained polyhedron V' = cone{w^ , . . . , v' n }, with v[ = (vi, 1) for 
i = 1, . . . , k and v[ = (v^ 0) for i = k + 1, . . .n, is referred to as polyhedral cone. 

By duality, the problem of converting a description by homogeneous inequalities 
into a description by extreme rays is equivalent to the opposite conversion problem. 
So for simplicity we assume from now on that V C ~R d is a polyhedral cone given by 
a minimal (non-redundant) set of generators (extreme rays). If 

V = cone{«i, . . . , v n } 

we say that V is generated by v%, . . . , v n G R d . 

We want to find a minimal set {fi, . . . , f m } C (M. d )* with 

V = {x G R d : fi(x) > Q,i = l,...,m}. 

By choosing a suitable projection, it is possible to reduce the problem further to 
the case where V is full-dimensional and does not contain any non-trivial linear 
subspaces. For example, if Vi,... ,v n G M. d span a /c-dimcnsional linear subspace, 
we may just choose (project onto) k independent coordinates. The appropriate 
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projection (i.e. equations of the linearity space) can be found efficiently via Gaussian 
elimination. In other words, without loss of generality we may assume the Vi span 
R d (V is full-dimensional) and the linear inequalities fi span (R d )* (V is pointed). 
A face of V is a set {x € V : f(x) — 0} where / is an element of 

V* = {f E {R d )* : f(x) > for all x e V}, 

the polyhedral cone dual to V . Note that (V*)* = V . The faces of a pointed 
polyhedral cone are themselves pointed polyhedral cones. We speak of a /c-face, if 
its dimension is k. The faces form a (combinatorial) lattice ordered by inclusion, the 
face lattice of V . The rank of a face in the lattice is given by its dimension. Each 
face is generated by a subset of the generators of V and therefore it is uniquely 
identified by some subset of {1, ...,n}. In particular, the 0-dimensional face is 
identified with the empty set and V itself with the full index set {1, . . . , n}. All 
other faces of V are identified with some strict, non-empty subset F C {1, . . . , n}. 

We write F < F' for two faces of V with F C F' and dim F = dim F' - 1. Two 
/c-faces of V are said to be adjacent, if they contain a common (k — l)-face and are 
contained in a common (k + l)-face. In particular, two extreme rays (1-faces) are 
adjacent, if they generate a common 2-face and two facets are adjacent, if they share 
a common (d— 2)-face (a ridge). By the properties of a lattice, for two faces F\ and 
i*2, not necessarily of the same dimension, there is always a unique largest common 
face F = F\ n F2 contained in them, and a unique smallest face F' , containing 
both. Any sub-lattice [F : F'] of the face lattice consisting of all faces containing F 
and contained in F' is known to be isomorphic to the face lattice of some pointed 
polyhedral cone of dimension dimi 7 " — dimF. Note that this is a special feature 
of face lattices of polyhedra, which is not true for general lattices. In particular 
the diamond property holds for face lattices: Every sub-lattice of rank 2 carries the 
combinatorics of a 2-dimensional polyhedral cone, namely a face F of rank k — 1, 
a face F' of rank k + 1 and two faces F\ and F 2 of rank fe. The rank 2 sub-lattices 
are in one-to-one correspondence to pairs of adjacent faces. 

A polyhedral complex A is a set of polyhedral cones (the cells of A) satisfying 
the following two properties: 

(a) If V € A then every face of V is also in A. 

(b) For all Vi,Vj € A, Vi H Vj is a face of both V l and V r 

The facets of a polyhedral cone V form a polyhedral complex, the boundary complex 
of V . Polyhedral complex A' is a subdivision of polyhedral complex A if every cell 
of A is the union of cells in A' and every cell of A' is contained in some cell of A. 
A subdivision is called a triangulation if every cell is a simplicial cone, i.e. is the 
conic hull of exactly d-extreme rays. By the homogeneous embedding of polytopes 
discussed above, we may equally discuss boundary complexes, subdivisions, and 
triangulations for polytopes. Two faces (cells) of a polytopal complex have empty 
intersection exactly when the corresponding cones intersect only at the origin. 

3. Polyhedral Symmetries 

3.1. Groups acting on polyhedra. In this paper we are especially inter- 
ested in the case where some non-trivial group acts on the given polyhedron V. 
The combinatorial automorphism group of T>, generated by n extreme rays, is the 
subgroup of Sym(n) of all permutations (acting on ray indices {1, ... ,n}) which 
preserve the whole face lattice of V . Thus the combinatorial automorphism group 
acts not only on the generating rays, but also on the set of facets (inequalities), and 
more generally on the sets of faces of a given dimension. Moreover, it preserves the 
inclusion relation among faces. It is known that the combinatorial automorphism 
group of V can be computed from the incidence relations between extreme rays and 
facets (see |KaSc03] ). 
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Given generators v\ , . . . , v n £ M. d of a polyhedral cone V and a subgroup G of 
the combinatorial automorphism group of V , we want to obtain the facets of V up 
to symmetry, that is, one representative F C {1, ... ,71} from each orbit under the 
action of G. 

If no group or only a small group is given, we usually still want to exploit as 
much symmetries of the given polyhedron as possible in the representation conver- 
sion. That is, ideally we would like to compute the full combinatorial automorphism 
group. However, we do not know how to compute it without computing the facets, 
which is precisely the problem we want to solve. So we have to settle for a com- 
promise and work with a more restricted type of automorphism. Note though, that 
one may obtain the full combinatorial automorphism group of V, after the facets 
have been computed 

In many cases the combinatorial automorphism group (or some nontrivial sub- 
group) may reflect geometric symmetries, for example if it has a representation as a 
subgroup of GLd(K) acting naturally on M. d and the polyhedral cone V . The group 
of all matrices in GLd(K) preserving V is called linear automorphism group of V. 
Since a linear automorphism permutes the set of extreme rays {R>o«i, ■ • • ,M> v n }, 
we naturally obtain a representation as a permutation group G < Sym(n). It is 
important to note here that although the induced permutation group is finite, the 
linear automorphism group of V is not necessarily so, and it can be quite awkward. 
Think for example of the quadrant in K 2 generated by the two non-negative coor- 
dinate axes. The induced permutation group is Sym(2), its linear automorphism 
group however is 

{(o J)'(d o) :a ' 6 < c ' deM > }- 

As explained in Section [2] we may limit our discussion to the special case of 
a full-dimensional, pointed polyhedral cone V generated by Vi,...,v n £ R d . The 
elements A £ GLd(R) of the linear automorphism group of V satisfy Av.i — AiU CT (i), 
where a £ Sym(n) is an induced permutation and > 0. Note, in case V is the 
homogenization of a (d — l)-dimensional polyhedron V' , the linear automorphism 
group of V corresponds to the so called projective automorphism group of V , that 
is, the set of all projective maps preserving V' . We are not aware of any practical 
algorithm to compute linear (or projective) automorphism groups or to decide if 
two polyhedral cones are linear (or projective) isomorphic. 

3.2. Restricted isomorphisms. A restricted isomorphism of two full-dimen- 
sional vector families V — {vi, . . . , v n } and V — {v[, . . . , v' n } in R d is given by 
a matrix A £ GLy(R) such that there exists a permutation a satisfying Avi — 
v 'a{i) f° r * = 1> • • • i n - A restricted automorphism of a vector family is a restricted 
isomorphism of V with itself. 

We speak of a restricted isomorphism between two polyhedral cones generated 
by the two vector families V and V , if it is a restricted isomorphism between V and 
V . Note though that the definition of restricted isomorphisms for polyhedral cones 
depends strongly on the choice of generators. Since such generators are only unique 
up to any positive factor, the choice of generators is very crucial. In practice the 
situation is usually not so bad as there is often a natural choice for the generators. 

The big advantage of restricted isomorphisms is that we can compute them by 
obtaining the graph isomorphisms for an edge colored graph. Given a vector family 
V = {vi, . . . , v n } C M d spanning M. d we define the positive definite matrix 

n 

(l) Q = £>4 

i—l 

Let the graph G(V) be the complete graph with vertices Vi and edge colours = 
vlQ~ 1 Vj. Then the following holds: 
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PROPOSITION 3.1. Let V,V C R d be two finite vector families. Then every 
isomorphism of the edge coloured graphs G(V) and G(V) yields a restricted iso- 
morphism of V and V' and vice versa. 

In practice the popular and very nice program nauty [MKa05 by McKay can 
be used to check for graph isomorphisms or compute the group of automorphisms. 
Note however, that the current version only takes vertex coloured graphs as input. 
We therefore have to transform our edge coloured graph to a somewhat larger 
vertex coloured graph that preserves the automorphism group (see for example 
[MKa06j . p.25). 

PROOF of 13. 11 Let Q be the matrix ([I} obtained from V. Denote by R the 
unique square root of Q , that is, the positive definite dx d matrix R with Q^ 1 = 
R 2 . Let Wi = Rvi for i = 1, . . . , n. Then the edge colours Cij of the graph G(V) 
are exactly the inner products w\wj of the transformed vectors W4 and Wj. In the 
same way we obtain Q', R' and w[ for V' . 

Now let A £ GLd(R) be a restricted isomorphism for V and V with associated 
permutation a £ Sym(n). Then we have AQA 1 = Q' . Moreover, the matrix 
T = R'AR^ 1 is orthogonal and satisfies Twi — ^(iV This implies = ^ a u) a (j)i 
i.e. the restricted isomorphism of the vector families V and V corresponds to an 
isomorphism between the edge coloured graphs G(V) and G(V). 

Suppose on the other hand a < Sym(rt) be an isomorphism between the edge 
coloured graphs G(V) and G(V'). By reordering elements, we may simply assume 
a = Id. Instead of looking for a solution of Avi — v[ we consider the equivalent 
equations Twi = w[. Since the Vi generate M. d , we find a basis (v i%1 . . . ,v id ) of R d . 
If P, respectively P' is the d x d matrix formed by (wi k ), respectively {w' ik ), then 
the matrix equation (cy) = (c^) takes the form P t P = P n P'. So, the matrix 
T = P'P^ 1 is orthogonal and we have for any k — 1, . . . , d and j = 1, . . . , n: 

w 'l k Tw j = {Twi^Twj = wj k Wj = w'l k w'y 

This yields w' i t (Twj — w' : j) — and since the form a basis of R d the relation Twj — 
w'j are implied. Thus we obtain a restricted isomorphism between V and V' . □ 

4. Orbits of faces 

4.1. Dealing with orbits. In order to generate facets (or more generally 
faces) of a polyhedron up to symmetries, it is necessary to deal with orbits of 
faces. In this section we briefly indicate what the basic tasks are that we have to 
accomplish and how these can be approached. 

As before, we assume that a polyhedral cone V is given by a set of generators 
vi,. . . ,v n £ M. d which define the extreme rays of V . Each face is represented by 
a subset of {1, ... ,n} which corresponds to the indices of generators incident to 
the face. We assume G < Sym(n) is some subgroup of the combinatorial automor- 
phism group of V, hence a permutation group acting not only on the set of indices 
{1, . . . , n} (respectively rays), but also on the whole face lattice of V . In particular, 
the dimension of a face and inclusion between faces are preserved by every group 
dement. 

In the most general form, the problem we have to solve is the following: Given 
two subgroups G\ and G2 of Sym(n) and a list L\ of Gi-inequivalent faces, we need 
to be able to obtain a list L 2 of G2-incquivalent faces. 

For example, if G\ = {Id} is trivial, the list L\ = {Fi, . . . , F^} simply is a 
list of pairwise unequal faces and we may obtain L2 by testing for G2-equivalcncc: 
Starting with L2 = {-Fi} and then subsequently adding Fi for i = 2, . . .,n to L2, 
if it is not G2-cquivalent to any element in L 2 . Clearly, this orbit fusion can be 
applied whenever Gi < G2. 
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In case G2 < G\ it is necessary to "break some symmetry" and split (factorize) 
orbits. This can be done with the double coset decomposition, described in Section 



4.2 Note, if neither G\ < G2 nor G2 < G\, one could in principle convert the G\- 
orbits into G2-orbits in two steps, by either using the intersection group G = GilHG^ 
or the group G = (Gi, G2) generated by G\ and G2 and obtaining G-orbits in an 
intermediate step. 

Orbit fusing and splitting is a typical and essential task when generating dis- 
crete struc tures up to isomorphism (see for example [GLM97 , [Ke99], [BrOOj . 



KaOs06 ) 



4.2. Fusing orbits, equivalence tests and canonical representatives. 

A common task is to decide whether or not two faces are G-equivalent, that is, 
whether or not the corresponding subsets of {1, . . . , n} lie in the same orbit under 
the action of G. 

The dimension of a face (rank in the face lattice) and its cardinality arc quickly 
testable invariants. Other invariants can easily be obtained, for example by look- 
ing at the action of G on pairs, triples, or other fc-tuples of indices (generators), 
respectively on lower dimensional faces. The number of elements from each such 
orbit included in a face is a G-invariant. Unfortunately, there is no clear rule how 
many such sets have to be chosen, so one has to rely on heuristics. 

If we only consider restricted isomorphisms, metric invariants can be used as 
well, for example the set of pairwise inner products between generators, discussed 
further in Section [7. 4. 21 

When all invariants are satisfied, then group computations must be done. For 
small groups it is possible to simply generate the whole orbit of a face. If the 
size of available memory is a problem, we may just keep a canonical representative 
for the orbit, for example the lexicographical minimum (when viewed as a sub- 
set of {1, . . . , n}). These can be found for example by a backtrack method (see 



Section 4.4 1 



4.3. Breaking symmetry by splitting orbits using double cosets. As- 
sume G2 < Gi < Sym(n) and that a list of Gi-orbits is given. In order to obtain 
a list of G2-orbits, we may split each orbit G\F, with a representative face F, by 
the well known double coset decomposition (cf. for example in [BrOO and [ Ke99p : 
The group G\ can be decomposed into double cosets 

r 

G 1 = \jG 2 g i Stab(G 1 ,F) 

i=l 

where gi, . . . ,g r are elements of G\. Then the orbit G±F is decomposed into 

r 

G\F = [J G 2gi F, 
i=l 

hence into r orbits G2-Fi with F, = giF. 

Note that the more straightforward algorithm of generating the full orbit or of 
computing a decomposition of G\ into cosets G2<?i is slower and/or requires more 
memory (see [Ke99]). 

4.4. Data structures and implementation issues. The fundamental data 
structures used to work with permutation groups are bases and strong generating 
sets (BSGS) (see for example [Se03 , HEO05J). Based on them a backtrack search 
on cosets of a point stabilizer chain can be used to obtain canonical representatives, 
to decide on (non-)equivalence as well as to obtain stab ilizers of subsets of {1, . . . , n} 



(faces). For details we refer to |Se03] and KaOs06 . An elaborate version is the 



partition backtrack introduced by Leon |Le91| (cf. [Le97 , [Se03j). These methods 
are known to work quite well in practice, although from a complexity point of view 
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the problems are known to be difficult. That is, there is no (worst-case) polynomial 
time algorithm known to solve these problems. Even worse, it is somewhat unlikely 
that there exist polynomial time algorithms, since the graph isomorphism problem, 
not known to be in P, is reducible to them in polynomial time (cf. Lu93 ). 

The computer algebra system GAP [GAP05 provides functions for generation 
of full orbits (Orbit) stabilizer computations (Stabilizer) and equivalence tests 
(RepresentativeAction) . 

In the case where we have special knowledge of the group G or its representation, 
it might be much easier to obtain canonical representatives or to compute stabi- 
lizers of faces. For example, we may have a situation where the symmetric group 
Sym(n) acts on n elements (see [An03 and [ChRo99j). Another example where 
the action of the symmetric group is used is described in |DeIn07] (see [Deln06 
for corresponding computer code). Typically, polyhedra arising in Combinatorial 
Optimization are convex hulls of (0/l)-vectors, where each coordinate (variable) 
represents an edge of a complete directed or undirected graph with n vertices on 
which Sym(n) acts. In [ChRe02 a method for obtaining canonical representatives 
in this situation is described. 



5. Decomposition methods 

In this section we describe two methods which reduce the facet generation 
problem under symmetries to a number of smaller instances of facet generation 
problems. In contrast to the original problem, solving the smaller problems (for sub- 
cones) may be feasible for available software, such as edd or Irs. These techniques 
have been proven to be successful in practice, in cases where standard approaches 
failed. 

We mainly distinguish two approaches, the Incidence Decomposition Method 



(see Section 5.1 1 and the Adjacency Decomposition Method (see Section 5.2 1. Both 



methods are reasonably natural and have been used separately by different authors. 

5.1. Incidence Decomposition Method. The Incidence Decomposition Me- 
thod reduces the problem of facet generation to a number of smaller problems, in 
which we generate facets that are incident to some extreme rays. As before let V 
be a polyhedral cone in K d , generated by Vi, . . . , v n . Assume G < Sym(n) is some 
permutation group acting on the face lattice of V . The set of extreme rays (indices) 
falls into orbits under the action of G. For each orbit we consider a representative 
Ti (index i) and generate a list of G-inequivalent facets of V incident to it. Then, 
in a post-processing step the lists of facets obtained in this way are merged to a 
list of G-inequivalent facets of V (see Section [4} . Since every G-orbit of facets of V 
contains a facet which is incident to one of the chosen representatives, the resulting 
list is complete. 

Input: n extreme rays of a polyhedral cone V and a group G < Sym(n) acting on 
"P's face lattice. 

Output: complete set T of G-incquivalcnt facets of V '. 
T ^- 0. 

1Z «— complete set of G-inequivalent extreme rays of V. 
for r e 1Z do 

T r «— facets of V incident to r. 
for F £ T r do 

if F is G-inequivalent to facets in T then 

T <- TU {F}. 
end if 
end for 
end for 
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The main computational gain comes from the following: when we compute 
the facets that are incident to a given ray fj, we may not have to consider all n 
extreme rays of V, because some of the rays may not be incident to facets which are 
incident to 7"j. These are exactly the extreme rays Tj = M.>QVj which give redundant 
inequalities of the polyhedron 

V* := {/ € (R d )* ■ f(vj) > 0,j = l,...,nand /(«<) = 0}. 

So, for each of the n — 1 rays with index in {1, . . . , n} \ {i} we may solve a linear 
program in (K rf )* to decide redundancy. As an outcome we obtain a list of n' < n 
rays needed in the definition of V* . The smaller n' , the bigger the computational 
gain. Note, that V* is of dimension d — 1. Its dual contains the line {Xvi : A G M.} 
and we may project it (or at least think of it as projected) down along this line to 
a d — 1 dimensional polyhedral cone. 

So the problem of enumerating facets incident to a given ray is a facet enumer- 
ation problem, but in one lower dimension and with fewer extreme rays. Some of 
the lower dimensional subproblems may still be too difficult; in this case we may 
apply the method recursively. We come back to this in Section [573) 

The Incidence Decomposition Method has been used in |Gr92] for finding 
the vertices of the metric poly tope MET 7. The method was also introduced in 
|ChRe9 6 along with the Adjacency Decomposition Method but it was found to be 
less competitive for their application. The Incidence Decomposition Method was 
also discussed by Fukuda and Prodon [PP96 . 

5.2. Adjacency Decomposition Method. As with the Incidence Decom- 
position Method, the Adjacency Decomposition Method is a reasonably natural 
method for computing the facets of a polytope up to symmetries. So it is no 
wonder that the method was discovered several times, for example, in [Ja93 as 
"algorithm de l'explorateur" , in ChRe96] as "adjacency decomposition method" 
and in DFPS01 as "subpolytope algorithm". Other example of applications are 
in |DuVa05j . [DeDu03] or |DDP02"] . The adjacency decomposition scheme is 
also implicit in polyhedral decomposition schemes, i.e. when a space is decomposed 
as an union of polyhedral cones. Examples are in |SY04j . [V0O8 and |Vo09| (cf. 
[DSV06bj and |DSV06ap . 

Rather than focusing on the incidence of (orbits of) facets to extreme rays as 
the Incidence Decomposition Method does, the Adjacency Decomposition Method 
focusses on the incidence of facets with other facets. Starting from a (set of) initial 
G-inequivalent facets, it traverses the adjacency graph of facet orbits. The initial 
facet(s) may be obtained by suitable linear programs, or in strongly polynomial 
time using the methods described in [BFM98 . 

Input: n extreme rays of a polyhedral cone V and a group G < Sym(n) acting on 
Vs face lattice. 

Output: complete set F of G-inequivalent facets of V. 
T <- {F} with F a facet of V. 
F <- 0. 

while there is an F e T do 
F ^FU{F}. 
T^T\{F}. 
H <— facets of F. 
for H e H do 

F' ^— facet of V adjacent to F along H . 

if F' is G-inequivalent to all facets in F U T then 
T^TU{F'}. 

end if 
end for 
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end while 

The facet F' of V with F fl F' = H for a given ridge H < F can be found 
by a gift-wrapping step (cf. |CK70j . |Swa85| ). Let vi,...,v n be the generators 
of V's extreme rays. The defining inequality / e (M d )* of the facet F' should 
satisfy f(vj) = for all generators u,; G F. The vector space of such functions has 
dimension 2. Let us select a basis {/i, /2} of it. If / = + 012/2 is the defining 
inequality of F, H or F', then f(vi) > for all i € {1, ... , n}. This translates into a 
set of linear inequalities on a.\, 02 defining a 2-dimensional pointed polyhedral cone. 
One finds easily its two generators (a\, a 2 )i<i<2- The corresponding inequalities 
fi = a \fi + "2/2 € define the two adjacent facets F and F' of V. In the 

special case where no d + 1 extreme rays lie on a hyperplane, the gift wrapping step 
corresponds to a simplex pivot (see |Chv83j ). 

A nice feature of the Adjacency Decomposition Method is that the adjacencies 
of the most symmetric facets, which are usually the most difficult to treat, may not 
have to be considered. This is due to the well known Balinksi's Theoreom (see e.g. 
[Zl98] . Theorem 3.14): 

Theorem 5.1. ^ |Ba61] ) Let V be a d- dimensional, pointed polyhedral cone. 
Let G be the undirected graph whose vertices are the facets of V and whose edges 
are the ridges ofV. Two vertices F±,F2 are connected by an edge E if F\C\F2 = E. 
Then, the graph G is id — \)-connected, i.e. removal of any d — 2 vertices leaves it 
connected. 

Using Balinski's Theorem, we know that if the number of facets in unfinished 
orbits (i.e. those whose neighbors are not known) is less than d — 1, then any 
G-inequivalent facet in the unfinished set must be adjacent to some already com- 
pletely treated facet. This simple criterion has proven |DS V06b] . |DSV07c] to be 
extremely useful in dealing with examples arising in the geometry of numbers. 

5.3. Recursion. Both, the Incidence Decomposition Method and the Ad- 
jacency Decomposition Method reduce the facet enumeration problem for a d- 
dimensional polyhedral cone V to a number of facet enumeration problems for 
cones in dimension d—1. These lower dimensional problems may be too difficult to 
treat with a standard method as well and therefore we might apply the Incidence or 
Adjacency Decomposition Method to these lower dimensional problems recursively. 
So we may speak of the Recursive Incidence Decomposition Method, the Recursive 
Adjacency Decomposition Method and the Recursive Decomposition Method, if a 
mixture of both is applied. 

The Recursive Incidence Decomposition Method has successfully been used for 
computing the vertices of the metric polytope METg in [DF MV03] , whereas the 
Recursive Adjacency Decomposition Method has successfully been used in |Ja93] 
and |DSV06bj . To the best of our knowledge, and to our surprise, a recursive 
mixture of both has not been used so far. 

The crucial steps for the recursion are "F r <— facets of V incident to r" (where 



r is some extreme ray) for the Incidence Method (see Section 5.1) and "F 



facets of F" (where F is some facet) for the Adjacency Decomposition Method 



(see Section 5.2 1. In both cases the problem is to obtain a list of facets for a (d— 1)- 
dimcnsional polyhedral cone F, whose extreme rays are given. Again we can exploit 
symmetries. As a result of a call of the Incidence or Adjacency Decomposition 
Method for F, we obtain a list of G^-inequivalent facets of F (ridges of V), where 
Gf is some group acting on the face lattice of F which we have to provide. In a 
post processing step we then have to obtain a list of G-p-inequivalcnt (respectively 
Stab(G-p, F)-inequivalent) facets of F (ridges of V) out of it. 

Note though, that the groups Gf and G-p may be unrelated, that is, the ele- 
ments of Cp do not have to be elements of G-p and vice versa. Think of examples 
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where V is a polyhedral cone without any symmetries, but with a facet having some 
symmetries. Vice versa, since we did not assume that Gf is the full (combinatorial) 
symmetry group, not even the stabilizer Stab(G-p, F) of F in G-p has to be a sub- 
group of Gf- In the latter case we may simply replace Gf by the group generated 
by Stab(G-p,F) and Gf, so that we assume Stab(G-p,F) < Gf- In this way we 
possibly enlarge the group and speed up the computations. Moreover, we are able 
to use the double coset decomposition (see Section 4.3) to split each G^-orbit GfF' 
of facets of F into a finite number of Stab(G-p, F)-orbits Stab(Gp, F)giF', where 
the gi € Gf represent the double cosets. 

Using the described decomposition methods recursively it might happen that we 
compute the facets of some sub-cones several times. For example if a face F satisfies 
F < Fx <\ V for a facet F\ of a polyhedral cone V, then there is exactly one other 
facet F2 such that F <3 F2 < V . Hence, if we apply the Adjacency Decomposition 
Method to F\ and F2, then we have to compute the dual description of F two 
times. Clearly, the number of such repetitions increases as the recursion depth 
increases. Moreover, equivalent sub-cones may occur in different parts of the face 
lattice. (Recall that any sub-lattice of the face lattice is the face lattice of some 
polyhedral cone.) To handle this, we propose to use a banking system. That is, given 
a difficult sub-cone V of which the facets have been computed up to symmetries, 
we store with V (represented by generators) its group of restricted automorphisms 
and a representative for each facet orbit. 

What about the recursion depth needed to practically solve a given polyhe- 
dral conversion problem? If we split the problem into subproblems with either 
decomposition method then some of the subproblems may be easy to solve and 
may not require a recursive treatment, others might be impossible to treat with- 
out such. The implementation |Dut07] allows one to provide a heuristic function 
to choose whether to recursively apply the Adjacency Decomposition Method, but 
in choosing the right level of recursion requires some trial and error. From our 
computational experience with this code (for example in [D uVa05j . |DeDu03] . 
[DDP02j . |DSV06aj . [DSV06b ), the incidence number of a face, that is, the 
number of extreme rays contained in it, gives a good measure for how difficult the 
representation conversion is for it. We therefore propose to treat the faces with low 
incidence numbers first and as much as possible without recursion. For subprob- 
lems to be potentially solved with pivoting methods (possibly up to symmetry as 
below) , the "probing feature" of Irs can be used to test for expected difficulties of 
a subproblem. 

Finally, let us remark, that parallelization is a possible way to speed up the 
(recursive) Decomposition Methods as well (see |ChRe02] and |DeIn06j ). 



6. An incremental method 

We briefly sketch here the Cascade algorithm by Jaquet [Ja93 , which is a sym- 
metry exploiting version of Fourier-Motzkin elimination. For this observe first that 
a d-dimensional polyhedral cone V generated by Vi,...,v n € M. d can be obtained 
as the projection of an n-dimensional polyhedral cone with n linear independent 
generators v[, . . . ,v' n . Without loss of generality, we may assume that (v±, . . . , Vd) 
is a basis of R d . We set v[ = ( Vl , n ~ d ) for i < d and v[ = (u i; Q"- d ) + e, for i > d. 
Let Pi(xi, . . . , x n ) — (xi, . . . , Xi, 0, . . . , 0) be the orthogonal projection of W 1 onto 
W x {0}™~\ Then for the cone V = cone-jX, ...,<} we have V = p d (T"). The 
cone V' is simplicial, since the family (f^)i<i<n is a basis of K". 

Take the rt-dimcnsional polyhedral cone V and the projection p n -i of K n on 
a (n — l)-dimensional hyperplane. The facets of the projection p n -i(V') are either 
projections of facets of V' or projections of intersections of facets of V '. In the 
Fourier-Motzkin elimination we first compute the facets of V and then successively 
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obtain the facets of p„_i (?-"), . . . , Pd('P') = V from one of these two cases. Note, 
since V is a simplex, the facets of V' are simply given by the (n — l)-subsets of 
extreme rays of V . 

In the Cascade algorithm we also consider symmetries in each step: If G is a 
group of symmetries of the polyhedral cone V, then the induced symmetry group 
on Pi(T") is the stabilizer of the set of vertices {vd+i, ■ ■ ■ , Vi] under G, which we 
denote by Gi. 

So, in order to compute the orbits of facets of PiiV') under d we need to 
compute the orbits of facets and ridges of pi^i(V'), first under Gi_i and then 
under Gi using the double coset method. 

It is well known that the Fourier-Motzkin method depends in a critical way on 
the ordering of generators v±, . . . , v n of V (see |ABS97] and |Bre99j ). The Cascade 
Algorithm introduces another such dependency as the size of the set-stabilizers vary 
enormously according to the chosen ordering. Among the many possible orders of 
the generators, it is not clear which ordering is best. 

In practice Fourier-Motzkin elimination suffers from generating many redun- 
dant inequalities. This problem can be eliminated by using the related Double 
Description Method [FP96 . For a description of how to use the double descrip- 
tion method to perform the Fourier-Motzkin projection steps without introducing 
redundant inequalities, see |Bre97| . Section 4.1. 

7. A pivoting method 

Since for certain classes of input, the most successful methods for the polyhedral 
representation transformation problem (without taking symmetries into account) 
are based on the pivot operation of the simplex method, it is natural to consider 
whether pivoting techniques can be adapted to the symmetric setting. 

7.1. The basis graph up to symmetry. For the purposes of this discussion, 
by k-basis, wc mean a set of k extreme rays generating the linear span of a fc-face 
of a pointed polyhedral cone. Two fc-bases are adjacent if they share k — 1 extreme 
rays. When k is not specified, we refer to [d— l)-bases of a d-dimensional polyhedral 
cone. Note that the bases of a polyhedral cone depend not only on the combinatorial 
information contained in the face lattice, but also on the linear dependencies among 
extreme rays incident to a given face. When discussing (d— l)-dimensional polytopes 
or their d-dimensional homogenizations, we use basis to refer interchangeably to a 
(d — l)-basis of the homogeneous cone or a (d — 2)-basis of the polytope (i.e., d— 1 
vertices spanning a facet of the polytope) . 

The basis graph has as nodes the bases, and as edges the pairs of adjacent 
bases. In the so-called non-degenerate case each facet of a d-dimcnsional cone is 
a [d — l)-dimensional simplicial cone. By a perturbation argument, for any cone 
there exists a triangulation of the boundary which is combinatorially equivalent to 
the boundary of a polyhedral cone with simplicial facets; it follows by Balinski's 
Theorem that the basis graph is at least (d — l)-connected. 

One of the earliest published methods of polyhedral representation transfor- 
mation [Cha53 is based on an exhaustive exploration of the basis graph. Using 
the pivot operation of the simplex method, adjacent (d — l)-bases can be found in 
time proportional to the number of input generators and the dimension (for details, 
see e.g. |AvOOj ). Pivoting methods based on Reverse Search [AF92 have mem- 
ory usage independent of the output size. In the present work we consider only 
more direct methods based on generation (and storage) of the basis graph up to 
symmetry. 

In the typical case, generating the entire basis graph is impractical, even if it 
is not stored. This is because of the enormous number of bases that correspond to 
each facet in the degenerate (i.e. non-simplicial facet) case. In the non-symmetric 
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setting, perturbation has been widely used to reduce the size of the basis graph 
under consideration, in particular via the lexicographic (symbolic) perturbation 
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discussed further in Section 7.3.2 We consider one way in which perturbation can 
be applied in the symmetric case in Section |7.3[ There is some tension between the 
notions of symmetry and perturbation though: whereas perturbation allows us to 
reduce the size of the basis graph, usually some of the given symmetries are lost. 
So there is a tradeoff in which ideally the quotient graph of the obtained new basis 
graph with respect to its remaining symmetries is as small as possible. 

7.1.1. Symmetry of the basis graph. The performance of pivoting based methods 
under symmetry is determined not by the total size of the basis graph, but by the 
number of orbits of bases. This number is not determined by the number of orbits 
of extreme rays and the number of orbits of facets. 

We define the basis automorphism group to be the subgroup of the combinatorial 
automorphism group that acts on the basis graph. To see that this is a non-trivial 
restriction, consider the cone generated by the following (row) vectors in M 5 



(2) 



Combinatorially, it is (the homogcnization of) a pyramid over an octahedron, and 
its combinatorial automorphism group is the 48 element octahedral group. The 
basis automorphism group has only 16 elements; in particular the orbit of the first 
d generators (as a basis) is of size 1. 

Example 7.1. The regular d-cross polytope Cd has one orbit of d-bases and one 
orbit of (d — l)-bases with respect to its automorphism group. 

PROOF. Since the d-cross polytope is simplicial, (i.e., its facets are all simplices) 
its (d — l)-bases are exactly its facets. Observe that at most one pair of vertices in 
a d-basis must consist of an opposite pair ±ej, since any two pairs are co-planar. It 
thus follows from the pigeonhole principle that a d-basis consists of a facet, along 
with one vertex not on that facet. All such simplices are equivalent under the 
automorphism group of the cross polytope. □ 

For a slightly more involved example, let V — conv(ui, . . . , v m ) C M. d and 
Q = conv(?Di, . . . , w n ) C R e . The wreath product of V with Q is defined as 

VlQ:= conv{(0,... ,0,^,0, ...,0,w k ) : 1 < i < m, 1 < k < n} C K nd+e . 

d{k-l) d(n-fc) 

So loosely speaking, the wreath product of V and Q is obtained by attaching 
to the vertices of Q pairwise orthogonal subspaces, each containing a copy of V . 

Example 7.2. The wreath product V I Q of a regular d-cross polytope V with 
a regular e-cross polytope Q has dimension D = 2de + e, 4de vertices, 2^ d+1 ^ e 
facets and one orbit of vertices, facets, and (D — \)-bases with respect to its linear 
automorphism group. 

Proof. The proof is based on the results of Section 2.2 in [ JL05] . The number 
of vertices of the cross polytopes V and Q are m — 2d and n — 2e. The assertion 
on the dimension and the number of vertices of V I Q follows immediately from 
the definition of wreath products. The number of facets and the number of bases 
orbits can be derived from Proposition 2.2 in |JL05| . By it, the facets of V I Q 
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are in one-to-one correspondence with choices (F; Fx, . . . , F e ) of a facet F of Q and 
facets Fi of V; which gives the count of facets. Assume w.l.o.g. that F contains the 
vertices v\, . . . , v e of Q. Then the rows of the following matrix describe the vertices 
contained in a facet of V I Q: 

V v 2 

V v e 
Fi v e+1 
F 2 v e+2 

F e V 2e _ 

The entries V and Fi stand for all choices of vertices from V, respectively Fj. So a 
facet of V I Q contains 3de vertices. 

The linear automorphism group of V I Q contains the semidirect product G = 
G P ><J Gq (where Gp and Gq denote the automorphism groups of V and Q and G P 
denotes an n-fold direct product). In particular, Gq is isomorphic to a subgroup 
of G, which permutes not only vi, . . . ,v 2e in the last e coordinates, but also the n 
copies of V attached to them. In contrast to this action interchanging the n d- 
dimensional subspaces, the elements of the subgroup G P of G act only within these 
subspaces. By this, all of the vertices and facets lie in one orbit under the action 
of G. The assertion on the bases follows as well, after noting that a base within a 
facet (as described above) is determined by e choices of (i-bases (in the e copies of 
V). All of them are equivalent with respect to G by Example |7.1| □ 



7.2. A prototype implementation. We have implemented a prototype called 
symbol ("Symmetry, Bases, and Lexicography") for facet enumeration up to sym- 
metries via pivoting. The code is mainly in GAP [GAP05 , with an external server 
based on Irs doing the pivoting. The main algorithm is a depth first search of 
the pivot graph; several options are provided for symmetry testing and perturba- 
tion (discussed further below). All symmetry groups are represented internally as 
permutation groups, and the backtracking algorithm of GAP's Representative- 
Action is used to test bases and facets for equivalence. Our main goal with this 
implementation is not to be able to attack large problems, but to provide a test-bed 
for better understanding of the main issues involved with pivoting under symmetry. 
In general we thus are not too concerned with a constant multiplicative slowdown 
in the runtime of our code versus more polished software; we are more interested in 
what kinds of problems can be solved in a reasonable amount of time. 

At least for certain special cases, our prototype is able to solve quite large 
problems. Figure [T] compares the run time of symbal with cdd and Irs on the 
wreath products of cross polytopes discussed in Example |7.2| Of course both of 
cdd and Irs compute all of the facets rather than orbits. The times reported here 
are on a 3GHz Pentium IV with 1G of memory. 

We call a cone (or polytope) basis- simplicial if it has as few basis orbits as 
facet orbits. Evidently every simplicial polytope is basis-simplicial. For perhaps 
less contrived examples of basis-simplicial polytopes than Example |7.2| we mention 
the Dirchlct-Voronoi-cells (DV-cells) of the root lattices D3, D4, and Eg (convex hull 
of their shortest non-zero vectors) and some related universally optimal spherical 
polytopes (see |BBC+06 ). Their boundaries consist of regular cross polytopes 
and regular simpliccs only and they thus have at most two orbits of basis with 
respect to their symmetry group (cf. Example 7.1). 
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Figure 1 . Experimental results for wreath products of cross polytopes 
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48 
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1440 
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4 
17 
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9892 
> 209000 



Table 1. Basis orbits of cubes 



As with the non-symmetric case, for cones with a large number of (orbits of) 
bases, pivoting is not a good approach, at least not without some kind of per- 
turbation. One benchmark of how orbitwise degenerate a polytope is, is whether 
the number of orbits of bases is larger than the number obtained by triangulat- 
ing the boundary. In this case, one is almost certainly better off applying either 
edd or Irs. A familiar example where this level of degeneracy occurs is in the 
(i-dimensional cubes (see Table [TJ . 

7.3. Orbitwise Perturbation. The basis automorphism group could in prin- 
ciple be larger than the linear automorphism group (as a simple example, consider 
a simplicial polytope combinatorially equivalent to the regular ri-cross polytope, 
but with trivial linear symmetry group). Nonetheless, for the reasons articulated 
in Section [3] and because linearity simplifies the discussion here, we consider here 
the case where we are given a subgroup of the restricted isomorphism group of a 
polyhedral cone. 

In this section we consider modifications of the standard lexicographic pertur- 
bation schemes that preserve some, but not necessarily all of the symmetry of the 
input. We first characterize one kind of transformation that preserves a prescribed 
symmetry group. 

Proposition 7.3. Let V be a vector family with restricted automorphism group 
G. Let H < G. Let V\, V2, ■ ■ ■ , Vk be the orbits of V under H . Let ui, . . . ,Uf. be a 
family of fixed points for H. Let 



V = \J{vi + Uj I Vl e V,} 



Let H' be the restricted automorphism group ofV. Then H < H' . 
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Proof. Let T £ H be a restricted automorphism for V and p the correspond- 
ing permutation such that Tv,i = i>p(j). Suppose vi £ Vj. Since there exists an 
automorphism mapping Vi to w p (i), by definition v p ^ £ Vj. It follows by linearity 
of T that 

T{v'i) = T(v t + Uj) = Tvi + Tuj = v p(l) + Uj = v' p{i) . 

□ 

Of course, in addition to preserving part of the symmetry group, we need to 
solve the original representation conversion problem. Let !/(•) denote the orbitwise 
perturbation map. The key property we need is that for each facet F of the original 
cone, there is some basis B such that v(B) is a basis of the perturbed cone. The 
standard way of ensuring this, and the one we adopt here, is to insist that the facets 
of the perturbed cone induce a subdivision of the facets of the original cone. Several 
notions of perturbation exist in the literature, including lexicographic |DOW55] , 
numeric MC89 . symbolic |Yap90 , and geometric [Se i98j . We will discuss a 
method that can be seen as a linear perturbation in the language of Seidel [Sei98 
or a modified version of the lexicographic perturbation first proposed by Dantzig, 
Orden and Wolf |DOW55j . 

Here we take a somewhat weaker definition of perturbation than is typical, 
since we are not concerned necessarily with obtaining a cone with simplicial facets 
as result. Let V be vector family. We say W C V is extreme (for V) if it is contained 
in some facet of cone{V). 

Definition 7.4. We say that V is a valid perturbation of V if there is a 
bijection v(-) between V and V such that for any W C V , 

(1) If u(W) is linearly dependent then W is. 

(2) If v{W) is extreme for V then W is extreme for V . 



From Definition 7.4 we have immediately the following refinement property. 



Proposition 7.5. For a valid perturbation V of a vector family V we have 

(1) The boundary complex of cone(V) is a subdivision of the boundary complex 
of cone(V) . 

(2) If X is a valid perturbation of V , then X is a valid perturbation of V . 

The following proposition shows that valid perturbations result from sufficiently 
small changes to a vector family. 

Proposition 7.6. For any vector family V C R d , for any W C V and any 
vector u £ R d , V(s) = { w + eu : w £ W} U V \ W is a valid perturbation of V for 
all e with |e| sufficiently small. 

Proof. Consider the vector family V{e) changing as e varies. From the point 
of view of linear dependence, we need only concern ourselves with a d-set of vectors 
v%, . . . , Vd, which are independent in V, but whose perturbation becomes dependent 
for e > 0. The following continuity argument shows that for small enough e no new 
linear dependencies may occur. Let T be the matrix whose rows are v\, . . . By 
renumbering we may assume the first k rows of T are in W . Let U be the matrix 
whose first k rows are u and the remaining rows are the zero vector. Then 

g(e) = det(T + eU) 

is a polynomial in e, which is non-zero in a neighborhood of e = 0. Taking the 
intersection of all such intervals over all full rank d-sets of vectors yields (1) and (2) 
of Definition WM □ 
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Let V be a vector family with restricted automorphism group G. Let H be 
a subgroup of G. Let Uh be a fixed point for H. We say that V is obtained by 
pushing W C V (respectively pulling W C V) if 

y- = { w + creM H |jueW}uy\W 

where a = 1 (respectively u = — 1) and e > is sufficiently small so that V is a 
valid perturbation of V. 

Let Vi, V2, ■ ■ ■ Vk be the orbits of V with respect to H. In general the com- 
binatorial structure of the resulting boundary complex depends not only whether 
each Vj is pulled or pushed, but on the order these operations are carried out. We 
say that V is an orbitwise lexicographic perturbation of V with respect to H if it is 
obtained by pulling or pushing each orbit defined by H in some perturbation order 
7T. Considering the special case where V contains the homogenization of the vertices 
of a d-polytope and uh is the homogenization of the origin, pulling (respectively 
pushing) corresponds to scaling an orbit outward (respectively inward) with respect 
to the origin; this terminology is consistent with that of Lee (Lee91 . 

Combining Propositions |7.3| |7.5[ and|7.6[ we have: 



Proposition 7.7. If V is an orbitwise lexicographic perturbation of V with 
respect to a subgroup H ofV's restricted automorphism group, then H acts on the 
basis graph of cone(V). 

In the case that there is only one orbit under the basis automorphism group, 
perturbing all input vectors by the same vector will not decrease the degeneracy of 
the problem. Our general strategy will thus be to choose a subgroup H that has 
multiple orbits of extreme rays, and at the same time choose different perturbations 
Ui for each orbit of V with respect to H. 

It is known [Lee91j that lexicographic perturbation of all of the vertices in- 
duces a triangulation (also called lexicographic) of the boundary of the polytope. 
In our case because we perturb (i.e. push or pull) all vertices in an orbit by the 
same amount, we cannot in general guarantee a simplicial result (i.e. an induced 
triangulation). Thus we continue to explore the complete basis graph of the per- 
turbed cone (polytope), up to equivalence classes of bases. We nonetheless hope 
for a significant reduction in the size of the basis graph by (effectively) breaking 
up very degenerate facets of the cone. The tradeoff, examined further in the next 
section, is that although we may lose some symmetry of the polyhedron, the quo- 
tient graph of the new basis graph (with respect to the remaining symmetries) may 
become smaller. Moreover, the degree of the vertices in the basis graph may have 
decreased, which speeds up the computation of the quotient. 

7.3.1. Choosing a subgroup to preserve. For any polyhedral cone V with sym- 
metry group G, in order to effectively use orbitwise perturbation, one needs to find 
a subgroup H < G such that the above mentioned possible computational gain is 
as large as possible. In order to develop some heuristics for how one might find 
such a subgroup, we have systematically studied the subgroups of the restricted 
automorphism group of the DV-ccll of the E7 root lattice (hereafter we use conv E7 
to denote this DV-cell) . Although we do not claim the results from a single example 
are in any way conclusive, they do at least suggest some ideas for further study. 

Let G denote the restricted automorphism group of conv E7. This is a group of 
order 2903040 that has one orbit on the 126 vertices of E7. The polytope conv E 7 
has 632 facets in 2 orbits. It has 161 basis orbits, compared to 20520 bases in the 
triangulation produced by Irs. 



The experiments on this section were carried out an Acenet (http://www. 
|ace-net . ca[ ) cluster with SunFire x4100 nodes (two 2.6 GHz dual-core Opteron 
285 SE processors and 4 GB RAM per core). 
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FIGURE 2. Number of basis orbits for various subgroups on conv E7 
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Figure 3. CPU time to enumerate facets for various subgroups on conv E7 



A sample of conjugacy classes of the subgroup lattice was generated using the 
GAP function LatticeByCyclicExtension, with the restriction that groups of size 
less than or equal to 100 were discarded. A representative was chosen from each 
conjugacy class, yielding 102 subgroups of G. In general combinatorially distinct 
perturbations could result from choosing a different perturbation order, and by 
varying the direction (push or pull) that each orbit is perturbed. In this experiment 
we restricted ourselves to orbitwise pulling, and ordered according to the smallest 
index in an orbit. Orbitwise pulling perturbations were generated for each of these 
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subgroups, and these perturbations were used as input to our prototype symbal. A 
total of 99 of these computations completed in the time allocated. The fastest time 
was about 16 seconds, and longest just over 26 hours. The unperturbed version 
completed in about half an hour. Perturbation thus yielded results ranging from 
more than a 100-fold speedup to more than 50-fold slow down. 

Figure [2] presents the experimental results in terms of the number of basis orbits 
computed. The groups are classified according to their order (rc-axis) and number 
of orbits of vertices ( shown by symbol). Recall that for any actual perturbation, 
we need at least two generator orbits. Of the subgroups with less than 64 orbits of 
(perturbed) bases all had two generator (i.e. vertex) orbits. The best result is 32 
basis orbits, achieved by a subgroup of order 80640. Within the class of subgroups 
having two generator orbits, the best results were achieved by groups of large order. 
Note that a large subgroup inducing two (or any other fixed small number) of 
generator orbits can be found by choosing random sets of elements as generators. 

Comparing Figure [2] with Figure [3] we see runtime is not completely a function 
of the number of basis orbits, but the size of the group also plays a role. The 
actual best time is achieved by a group of order 336, with 4 generator orbits and 
82 basis orbits. After a certain point a small subgroup is no longer advantageous, 
since the number of basis orbits are too large. In the limiting case, with a trivial 
group, our pivoting scheme produces a lexicographic triangulation like Irs; albeit 
less efficiently in our prototype implementation. A relatively small subgroup with 
a fixed number of input orbits can also be found by the same random sampling 
strategy. 

Figure [4] illustrates our two heuristics applied to the three dimensional cube. 
On the left we have the triangulation induced by a subgroup of order 24 having two 
generator orbits, In the center, the triangulation is induced by the stabilizer of a 
pair of opposite vertices. On the right, we take the same stabilizer, but push the 
first orbit instead of pulling. 

The triangulation induced by pulling opposite vertices of a 3-cube turns out to 
be a special case of a triangulation of the d-cube that has only one basis orbit in any 
dimension. Let I d = [— 1,1] denote the centrally symmetric d-cube. Define e = 
Si=i e i For each permutation p € Sym(d), let A p denote the simplex with vertices: 

-e, -e + 2(e p( i)), -e + 2(e p(1) + e p(2 )), . . . , e 

Let A denote the union of all A p and let A denote the set of (d— l)-simplices formed 
by intersecting the simplices of A with the boundary of I d . It is known [DRS07 
that A forms a triangulation of I d ; consequently A forms a triangulation of the 
boundary, since every triangulation of a polytope induces a triangulation of its 
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boundary. We call A (respectively A) the linear ordering triangulation of I d (re- 
spectively of the boundary of I .) Since there is a bijection between permutations in 
Sym(d) and simplices in A, Sym(d) acts transitively on A by permuting coordinates. 

Example 7.8. Let V denote the vertices of I d , i.e. V = {±l} d . Let H de- 
note the stabilizer of the automorphism group of L d on {— e,e}. Define uj(v) = 
mm(e T v , —e T v) . Let v(V) denote the H-orbitwise pulling of V in the order in- 
duced by uj. Then the following holds: 

(a) H acts transitively on A. 

(b) dcomfv{V) is combinatorially equivalent to A. 

Proof. Let us first remark that the order of orbits induced by 10 is well defined. 
The group H < Sym(2 d ) is generated by a set of generators permuting coordinates, 
along with a switching permutation a that maps v to —v. It follows that the H- 
orbits of V are the equivalence classes of V with respect to uj. 

To see (a), consider the d-simplex (corresponding to the identity permutation) 
Aid with vertices { — e, e + 2e\, — e + 2(e\ + e2), . . . ,e}. Let a' denote the per- 
mutation in H that first applies a (i.e. switching) followed by reversing the order 
of coordinates. The permutation a' is an automorphism of Aid which carries the 
(d — l)-simplex { e, — e + 2e l7 . . . , e — 2e^) } to { e,e — 2e^, . . . , — e + 2e 1 }. These 
two (d — 1) simplices are precisely the contribution of Aid to the linear ordering 
triangulation A. Thus any simplex of A can be mapped to any other by an action 
of Sym(d) on A (i.e. permuting coordinates), followed by possibly applying a' . 

We now consider (b). We argue that v(V) induces a linear ordering triangula- 
tion of each fc-face of L d , < k < d. Each (d — 2)-face will receive the same (linear 
ordering) triangulation from the two facets that contain it, hence the triangulations 
of the facets form a triangulation of the boundary. 

For < k < 2, there is nothing to prove. Let F be a fc-face of L d , 2 < k < d. Let 
v + (resp. v~) be the vertex of F with the most positive (resp. negative) coordinates. 
Recall that we will first pull the ii-orbit with smallest u> value. 

If the functional uj is minimized uniquely at v* G {v~,v + } then the pertur- 
bation corresponds locally to a standard [Lee91J pulling and the corresponding 
subdivision is into pyramids F±, . . . ,Fj with apex v* and bases corresponding to all 
of the k — 1 faces of F that do not contain v* . 

Otherwise uj is minimized at both v~ and v + . The perturbation thus takes v + 
to pv + and v~ to some pv~ , p > 1. This turns out to induce a subdivision of F into 
poly topes F? with vertices Vj — { pv~ , pv + } U Rj where Rj is the vertex set of a 
(k — 2)-face of F containing neither v~ nor v + . The polytope F? is a 2-fold pyramid, 
since v~ ^ aff Rj and v + aff({ v~ } U Rj). It follows that dimFj = dimF = k. 

That the F? are induced by the perturbation can be seen by exhibiting a sup- 
porting hyperplane of conv v(V). Without loss of generality, let F be defined by 
equations Xi = 1, i — 1, . . . ,d — k. The (fc — 2)-face conv(Rj) must be defined 
by further equations x p = 1, x q = —1, p, q > d — k. Let /i = 1/p. Consider the 
hyperplane hj = { x \ a T x — 1 }, where 

'fj,/(d-k) \<i<d-k 

(l-/x)/2 i=P 

-(l-/i)/2 i = q 

^ otherwise . 

It can be verified that hj supports conv v(V) and hj n conv^(y) = Ff. 

It remains to see that the Fj cover F, i.e. that there are no other cells in the 
induced subdivision. Consider an arbitrary relative interior point x of F. Let r be 
the ray from pv~ through x. Let y be the first intersection of r with dF after x. 
For each (k — l)-face of F, the double pulling of { v + , v~ } acts like a single pulling 
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decomposing the boundary of F into pyramids with apex either pv~ or pv + ; hence 
y € dFj for some j. It follows that x is in Fj. 

Now suppose for all j < k, v(V) induces a linear orderi ng t riangulation of 



the j'-faces of I d . From the refinement property Proposition 7.5 we know that 
v(V) induces a decomposition of the pyramids F±, . . .,Fk (respectively of the 2-fold 
pyramids Ff,... ,F^ k _^). In both cases the resulting /c-simplices correspond to 

the coordinatewise-monotone paths from v~ to v + in F. □ 

7.3.2. Symbolic implementation. Suppose we are given V' C R d+1 which are 
homogeneous coordinates for the generators V of some polyhedron V C M. d . We 
consider V here as a matrix (with Vi as rows) , and let u denote the column vector of 
corresponding (d+ l)st coordinates in V' . Following the conventions of Section [2] 
we suppose Uj is 1 if Vj is a vertex, and if it is an extreme ray. We consider the 
polyhedron V® 

(3) V° = { x £ R d : Vx > -u } . 

The polyhedron V may be thought of as cone(V)* C (M <i+1 )* intersected with 
the hypcrplane x^+i = 1. By duality, to find the generators of V is equivalent to 
finding the facets of V ■ 

Let G be the restricted automorphism group of V. Let G be the induced 
group acting on V . We will assume without loss of generality that the origin is the 
ccntroid of V and thus a fixed point of G. It follows that e^+i is a fixed point of 
G' . Applying Proposition 7.3 to our dual representation ([3]), we see that orbitwise 



perturbing the right-hand side vector according H < G will preserve the symmetries 
of H. The perturbed system thus has the form 

(4) Vx>-(u+p), 

where pi = o~jEj for <jj € {±1} and j is the index of the orbit containing Uj. To 
ensure an orbitwise lexicographic perturbation, we will insist 

(5) 1 > e x > e 2 > • • • > £k > 0, 

where by x 3> y we mean that y is much smaller than x, i.e. it is not possible to 
combinatorially change the polyhedron defined by Q by choosing y > smaller. To 
implement this symbolically, we need a modification of the standard lexicographic 
pivot rule (see [AvOOj or |Chv83j for more details). Let b = —(u + p) and A — 
[V —I]. After adding slack variables to Q, we are left with a system of the form 

Ax = b 

with n rows and n + d columns, where the first d columns are the decision variables. 

A feasible basis consists of a partition (/?, rf) of the column indices such that Ap 
(columns of A indexed by j3) is non-singular and any slack variables in x*^ = A^ 1 b 
are non-negative. In order to move from one feasible basis to another, we need 
to perform a pivot. We start by choosing a column index j to leave (3. To find a 
column index to replace j, we need to find 

. K 

argmuij — 

where b* = A~Z b and a = A^ 1 A v . In our case, where b = —u—p, we may decompose 
b* , and thus the ratio test into two parts A^" 1 (— u) and A^" 1 (— p). Because the values 
of the £j are chosen very small, the second part is considered only to break ties. We 
write 

£l 



V(-<") = N 



t2 
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Figure 5. 3 facets with a 3- fold rotational symmetry. Basis has 
vertices Pi,qi, around the boundary of the facet. 
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Figure 6. Comparing the search with pruning (right) and without (left). 



where column j of N is defined by summing the columns of Ap 1 corresponding to 
orbit j of generators, and multiplying by — <jj. Because of the ordering (j5j, in order 
to evaluate 



we proceed column by column in N, reducing the set of ties at each iteration. 

7.4. Other refinements and implementation details. 

7.4.1. Adjacency Decomposition Pruning. Consider facets Fq and Fi that are 
equivalent under some symmetry of the basis automorphism group. This same sym- 
metry acts as an isomorphism between the corresponding basis graphs. It follows 
that when we discover a basis B defining a new orbit, but the facet F spanned by B 
is known, we do not need to explore the neighbours of B since they will be explored 
in our canonical (i.e. discovered first) facet in the orbit of F. In order to ensure 
orbits are not discarded, we are careful not to mark B as known until its canonical 
discovery. 

Although this pruning does not reduce the number of orbits of bases explored, 
it can reduce the number of actual bases visited (and tested for isomorphism) , since 
bases of a given orbit are not revisited in every copy of the facet F. As an example, 
consider the 3 quadrilateral facets illustrated in Figure [5] with rotational symmetry 
yielding 4 orbits of bases. Without pruning, a depth first search visits all of the 
bases; with pruning only 8 of 12 bases are visited (see Figure ^J. The speedup 
obtained depends roughly on the number of facets visited by the unpruned search, 
which is bounded by the number of basis orbits. 

7.4.2. Metric Invariants. In the case where our symmetry group preserves the 
inner product between pairs of vectors, as is the case for the restricted automor- 
phisms discussed in this paper, we may take advantage of this in several ways. 
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For any face or basis X to be tested for isomorphism, we may construct a graph 
(analogous to that constructed in Proposition 3.1l whose nodes are the vectors of 
X and whose edges are the angles between them. This graph contains geometric 
information not present in the index sets representing X, which can help to speed 
up an algorithm to find an isomorphism. A simpler observation, and equally widely 
applicable, is that the set of pairwise inner products of two isomorphic faces or bases 
must be equal. This allows us to store orbit representatives in a data structure such 
as a hash table or a balanced tree, with the key to the data structure being the set 
of inner products. This permits more efficient isometry testing by retrieving exactly 
those orbit representatives which pass the inner product invariant. 

It is computationally easy to test whether a given linear transformation T is 
in the restricted automorphism group Aut(V) of vector family V . Since we are 
interested in restricted automorphisms carrying basis X to basis Y, we can addi- 
tionally test if HXT = Y for some II € Aut(A) (where Aut(A) can be computed by 
the same techniques as Proposition 3.1 ). We have only implemented an exhaustive 
search of Aut(A), and this is naturally only effective when Aut(A) is quite small. 
In principle it should be possible to integrate the test for T being a restricted iso- 
morphism into a backtracking procedure to search for II. 



8. Conclusions 

Much as in the case of polyhedral representation conversion without symme- 
tries, a certain amount of trial and error seems to be necessary to decide on the 
the best method to attack a given conversion problem up to symmetries. Currently 
decomposition methods have the best record of solving interesting problems; on 
the other hand current software requires a certain amount of user intervention in 
the form of choosing how to treat subproblems. It would be helpful to automate 
this process. In this context, a virtue of the pivoting methods is that good meth- 
ods to estimate their running time exist [AD94 . It would be beneficial, not just 
when working with symmetry, to have effective methods (or at least heuristics) for 
estimating the running time of incremental methods. 
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